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Abstract 



An ab initio method for calculating electron-phonon coupling parameters 
is presented. The method is an extension of the plane-wave-based linear- 
response method for the calculation of lattice dynamics. Results for the mass 
enhancement parameter A and the electron-phonon spectral function a^F(u;) 
for Al, Pb and Li are presented. Comparisons are made to available experi- 
mental data. 
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I. INTRODUCTION 



The electron-phonon interaction in metals plays an important role in a variety of exper- 
imentally accessible quantities including the enhancement of the electron mass, the phonon 
lifetime arising from electron-phonon scattering, electrical and thermal conductivities, and 
the superconducting transition temperature. The electron-phonon spectral function a'^F{(jj) 
measures the effectiveness of phonons with energy uj to scatter electrons from one part of the 
Fermi surface to another part. Once it and the Coulomb pseudopotential /x* are determined 
for a superconducting material, all of the thermodynamic properties of the superconductor, 
including the superconducting g function of temperature, the transition temperature, 

and the discontinuity in the specific heat at Tc, can be computed In addition, transport 
properties of materials in the normal state can be calculated from the closely related spectral 
function a%F{ijj). 

The ability to accurately calculate electron-phonon coupling parameters has long been a 
sought after goal. It is a formidable task requiring knowledge of the the low-energy electronic 
excitation spectrum, the complete vibrational spectrum, and the self-consistent response 
of the electronic system to lattice vibrations. Ah initio calculations of electron-phonon 
coupling parameters have proceeded along two distinct lines. In the rigid ion (RI) and 
rigid muffin tin (RMT) schemes, the screened electron-phonon potential is approximated 
by neglecting changes in the potential everywhere except within the atomic sphere of the 
displaced atom. While these non-self-consistent approximations appear to be adequate for 
many transition metals @], their validity has been questioned in some cases, especially 
for anisotropic or low-density-of-states materials 1^-0]. An alternative to the RI and RMT 
methods is the frozen-phonon total-energy method p|. In this approach, the electron-phonon 
matrix elements are evaluated using the self-consistently screened potentials corresponding 
to frozen-in phonon displacements. The primary drawback of the frozen-phonon approach is 
that only phonon wavevectors that are commensurate with the lattice and that correspond to 
reasonably sized supercells can be considered. This makes it difficult to determine accurately 
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quantities that involve integrations over the phonon wavevector q throughout the Brillouin 
zone. These include, for example, the electron-phonon mass enhancement parameter A, the 
phonon density of states F{uj), and the electron-phonon spectral function a'^F{uj). 

Recently, linear-response theory within the framework of density-functional calculations 
has been shown to be an efficient and powerful alternative to the frozen-phonon method for 
calculating lattice dynamical properties of solids §]. Atomic displacements are treated as 
perturbations, and the electronic response to the perturbation is calculated self-consistently. 
Perturbations of arbitrary wavevector q can be treated without using supercells. The linear- 
response method has been implemented with a variety of different basis sets for representing 
the electronic wavef unctions, and it has been successfully applied to the study of lattice 
dynamics in a wide range of systems including semiconductors PJTO[|, metals [|lT]-0, ferro- 
electrics |]14[ , and surfaces . 



In this paper, we extend the plane-wave-based density-functional linear-response method 
to the calculation of electron-phonon coupling parameters. A similar method based on 
linearized muffin-tin orbitals was recently presented in Ref. 16. The use of a plane-wave 
basis has the advantage of simplicity both in terms of formalism and implementation. We 
apply our method to the calculation of the electron-phonon spectral function a'^F{uj) for 
several elemental metals: Al, which is a well-studied weak-coupling superconductor; Pb, 
which is the prototypical strong-coupling superconductor; and Li, which is not observed to 
have a superconducting phase, but which earlier calculations have found to have a moderate 
electron-phonon mass enhancement parameter. 

II. METHOD 

For notational simplicity, we consider the case of a single atom of mass M per unit cell. 
The electron-phonon matrix element for scattering of an electron from a Bloch state nk to 
another Bloch state n'k.' by a phonon of frequency u^k-kv is 

g{nk,n'k',iy) = {-—^ )^^^(V^nk|ek-kV ■ VaKd^n'k'); (1) 



where ek-kv is the phonon polarization vector, and VrKc is the gradient of the self- 
consistent potential with respect to atomic displacements. The linewidth of phonon qt/ 
arising from electron-phonon scattering is given by 

7q^ = 2Txuj^y ^ 5(E„k - EF)6{En'\<i' - -Ei;^)5k-k'-q|fi'(nk, n'k', z/)p, (2) 

n,n',k,k' 

where ii^i;' is the Fermi energy. 

For many applications, the quantity of interest involves a sum or average of the electron- 
phonon coupling strength over wavevectors throughout the Brillouin zone. For example, the 
electron-phonon spectral function a'^F{ijj), which is a central quantity in the strong-coupling 
theory of superconductivity, is given by a sum over contributions to the coupling from each 
phonon mode: 

Here N{Ep) is the electronic density of states at the Fermi level. The dimensionless electron- 
phonon mass enhancement parameter also involves a sum over modes and can be expressed 
as the first inverse frequency moment of the spectral function: 

A = duja'^F{uj)/uj. (4) 

In this work, the electronic wavefunctions ipnk and eigenvalues Enk are calculated us- 
ing the ab initio pseudopotential local-density formalism. The electron-ion interaction is 
represented by soft separable pseudopotentials and the single-particle wavefunctions 
are expanded in a plane- wave basis set. The Wigner form of the exchange and correlation 
functional is employed |jl8|, and in the case of Pb and Li, the partial core correction is used 



to handle the nonlinearity of the exchange and correlation interaction between the core and 
valence charge densities |19|. Unless otherwise indicated, the calculations are performed 
at lattice constants that are determined within the local-density approximation (LDA) and 
that are in good agreement with the experimental values. 

The phonon frequencies and polarization vectors are calculated using linear-response 
theory. The second-order change in the total energy, and hence the dynamical matrix. 



depends only on the first-order change in the electronic charge density. The linear response 
of the electronic density to atomic displacements is determined self-consistently by solving a 
Bethe-Salpeter equation as discussed in Ref. 11. We have generalized the method to include 



corrections for the overlap between core and valence charge densities |]T2[. The electron- 
phonon matrix elements, g{nh,n'k.' ,1^), are easily computed from the first-order change in 
the self-consistent potential. 

The doubly-constrained Fermi surface sums in Eq. (2) are performed using dense meshes 
of 1300 and 728 k points in the irreducible Brillouin zones (IBZ) of the fee and bcc structures, 
respectively. The 6 functions in energy are replaced by Gaussians of width 0.02 Ry. Because 
of the large number of k points sampled, the results are not very sensitive to the Gaussian 
width. Phonon wavevectors are sampled on coarser meshes of 89 and 140 points in the fee 
and bcc IBZs, respectively. 

III. RESULTS 

To test the accuracy of the method, we consider first the simple metal Al. The LDA gives 
very good structural properties for fee Al, and the linear-response method yields phonon dis- 



persion curves in excellent agreement with experiments throughout the Brillouin zone [|lT 
The a'^F{uj) for Al calculated in this work is shown in Fig. 1 (solid line), along with re- 
sults from two experiments (long- and short-dashed lines) pOf. Extraction of a'^F from 
conventional McMillan-Rowell tunneling spectroscopy is not possible for Al since it is too 
close to an ideal BCS superconductor. Instead the experimental spectral functions shown in 
Fig. 1 were extracted from proximity electron tunneling data . Unfortunately, the inver- 



sion of this type of tunneling data involves the introduction of additional fitting parameters 
characterizing the proximity layer. This introduces uncertainties in the extracted spectral 
functions, as evidenced by the differences between the two experimentally determined a'^F 
functions shown in Fig. 1. The calculated spectral function agrees reasonably well with the 
experimental curves. 
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The value of the electron-phonon mass enhancement parameter determined from the 
first inverse-frequency moment of the calculated spectral function is A = 0.438, which is 
in good agreement with other linear- response and frozen-phonon calculations and 



with heat capacity data Within Eliashberg theory, Tc is a functional of the spectral 
function a'^F{uj) and the Coulomb parameter fi*. Using the calculated a'^F as input into the 
Eliashberg equations, we find that a /i* of 0.162 is needed to obtain the measured transition 
temperature of Tc=1.18 K. The same value of /i* yields a gap equal to the experimental 
value of 0.180 meV. This consistency between independent fits to the gap and to Tc is a 
measure of the accuracy of our results for a'^F. 

We consider next the case of the strong-coupling superconductor Pb, for which high- 
quality conventional tunneling data are available. From the theoretical stand point, the 
importance of relativistic effects in Pb make it a more difficult system to treat than Al. 
The present calculations for Pb are performed in the scalar relativistic approximation. The 
phonon dispersion curves are shown in Fig. 2. Overall, there is good agreement between 
calculated (solid lines) and measured |^ (circles) phonon frequencies, and the calculation 
is able to reproduce some subtle features in the dispersion curves such as the Kohn anomaly 
in the longitudinal mode along the F to K direction. Note however that there are significant 
quantitative discrepancies between the calculated and measured frequencies for the low- 
energy transverse mode, especially in the regions near X and K. The minima in both the 
longitudinal and transverse modes at X are not observed in other fee metals, and they 
suggest the presence of very-long-range forces. Indeed, as shown by the dashed curves in 



Fig. 2, an eighth-neighbor Born-von Karman fit |2J] to the measured frequencies is unable 
to reproduce the dispersion near X, especially in the case of the longitudinal mode. Over 
the years, a number of exotic mechanisms have been proposed to explain the unusual shape 
of the dispersion curves in Pb Our preliminary results obtained using the frozen- 

phonon approach indicate that differences between the calculated and measured frequencies 
are significantly reduced if the spin-orbit interaction is taken into account. 

The Eliashberg function for Pb is plotted in Fig. 3. The calculations (solid line) yield the 



two-peak structure seen in the data (circles) p6| , but there are differences in the peak loca- 
tions and heights, especially in the case of the lower-frequency peak. The mass enhancement 
parameter is calculated to be A = 1.20, which is significantly lower than the tunneling result 
of 1.55. These discrepancies are due in part to the errors in the calculated transverse-mode 
phonon frequencies. Note that the location of peaks in the spectral function is determined 
to a large extent by the location of peaks in the phonon density of states F{uj) {i.e., by 
the delta functions in frequency in Eq. (3)). The overestimation of the transverse-mode 
frequencies in our calculation results in an upward shift of the lower-frequency peak in both 
F{uj) and a'^F{uj). We have also computed the spectral function using the phonon frequen- 
cies generated from the Born- von Karman fit to the data along with the calculated phonon 
linewidths The resulting a'^F{uj) is plotted as a dashed line in Fig. 3. Using the 



empirical force constants, which accurately describe the dispersion of the transverse modes, 
we obtain good agreement with the experimental results in the low-frequency regime. On 
the other hand, since the Born-von Karman fit does not yield accurate frequencies for the 
longitudinal mode, the resulting spectral function is less accurate than the first-principles 
result in the high-frequency regime. It appears that in order to improve our description 
of the electron-phonon coupling in Pb, it will be necessary to modify the computational 
method to take into account the relativistic spin-orbit coupling interaction. 

Finally, we examine the electron-phonon coupling in bcc Li. The lack of a superconduct- 
ing transition in Li has been a long-standing puzzle. Both frozen-phonon p8[ and RMT ||2^ 



calculations have suggested that the electron-phonon coupling strength in Li is similar to 
that in Al. This would suggest a transition temperature on the order of 1 K if a value of 
/i* ^ 0.15 is assumed. Experimentally, however, no transition is observed, at least down to 
6 mK g. 

The LDA tends to underestimate the lattice constant of alkali metals. For Li, the cal- 
culated lattice constant of 3.41 A is about 2.3% smaller than the measured value of 3.49 
A. We have performed calculations at both values of the lattice constant. The phonon fre- 
quencies calculated at the experimental lattice constant are plotted in Fig. 4. Overall, the 



frequencies agree well with the neutron diffraction data |^ , and subtle features such as the 
crossing of the longitudinal and transverse modes along the F to H direction are reproduced 
by our calculations. At the LDA lattice constant, the calculated phonon frequencies increase 
by up to 8%. This is consistent with recent supercell calculations carried out at the LDA 
lattice constant determined without including core corrections (3.35 A) [^. In that case, a 
constant scale factor of about 0.86 was needed to bring the theoretical results in line with 
the experimental frequencies. 

The mass enhancement parameter in Li is calculated to be A = 0.45 and 0.51 at the exper- 
imental and (core-corrected) LDA lattice constants, respectively. These values are similar to 
results obtained earlier within the RMT approximation or using the frozen-phonon method. 
Using our result for a'^F{uj) as input into the Eliashberg equations, we have calculated the 
superconducting transition temperature as a function of the Coulomb pseudopotential. As 
shown in Fig. 5, an unphysically large value of /i* ~ 0.28 is required in order to suppress 
Tc below the experimental limit. The observed absence of superconductivity in Li therefore 
remains an open problem. The resolution of this puzzle may require consideration of the 
low-temperature crystal structure of Li, in which the electron-phonon interaction may be 
weaker PSI, as well as the role of manybody interactions such as electronic correlations and 



spin fluctuations p3 |. 



In summary, we have presented an accurate and efficient method for calculating electron- 
phonon coupling parameters from first principles. The electronic response to the atomic 
displacements is determined self-consistently and phonons of arbitrary wavevector can be 
treated. This method, which is an extension of the plane-wave-based linear-response method 
for calculating lattice dynamics, is applicable to a wide range of materials. Results for the 
mass enhancement parameter A and the electron-phonon spectral function a'^F{uj) for Al 
are in excellent agreement with available experimental data. In the case of Pb, there are 
larger discrepancies between our results and the tunneling data, but we attribute this to 
the neglect of the spin-orbit interaction in our calculations. The present results for the 
mass enhancement factor in bcc Li are in accord with earlier calculations, indicating that 
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further theoretical investigations are needed to resolve the question of the absence of a 
superconducting transition in this material. 
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FIGURES 

FIG. 1. Electron-phonon spectral function for Al. Results from the present calculation are 

represented by the solid line. Results from proximity electron tunneling spectroscopy experiments 
are indicated by the long- and short-dashed lines. 

FIG. 2. Phonon dispersion curves for Pb. The solid lines connect frequencies calculated at 
the sampled wavevectors, and the circles indicate the experimentally measured frequencies. The 
dashed line is an eighth-neighbor Born-von Karman fit to the measured frequencies. 

FIG. 3. Electron-phonon spectral function for Pb. The results based on the calculated frequen- 
cies are shown as a solid line, those based on the frequencies obtained from the force-constant fit 
are given by the dashed line, and the results from the inversion of tunneling data are plotted as 
circles. 

FIG. 4. Phonon dispersion curves for Li. The solid lines connect frequencies calculated at the 
sampled wavevectors, and the circles denote the experimentally measured frequencies. 

FIG. 5. Transition temperature vs. Coulomb pseudopotential for Li. The transition tempera- 
ture, which is plotted on a log scale, is computed as a function of /x* using the calculated a^F as 
input to the Eliashberg equations. 
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